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. Abstract 

\ Social, technological and economic time series are divided by events which are usually assumed 

f ' to be random albeit with some hierarchical structure. It is well known that the interevent statistics 

^ O^ . observed in these contexts differs from the Poissonian profile by being long-tailed distributed with 

! resting and active periods interwoven. Understanding mechanisms generating consistent statistics 

> 

Qf^ ' have therefore become a central issue. The approach we present is taken from the Continuous 

in 

en ' Time Random Walk formalism and represents an analytical alternative to models of non-trivial 
^ ■ 

ly-^ , priority that have been recently proposed. Our analysis also goes one step further by looking at 
O , 

00 , the multifractal structure of the interevent times of human decisions. We here analyze the inter- 
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transaction time intervals of several financial markets. We observe that empirical data describes a 



^ . subtle multifractal behavior. Our model explains this structure by taking the pausing-time density 

H : 

, 5t , in the form of a superstatistics where the integral kernel quantifies the heterogeneous nature of the 

executed tasks. An stretched exponential kernel provides a multifractal profile valid for a certain 
limited range. A suggested heuristic analytical profile is capable of covering a broader region. 
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I. INTRODUCTION 



The dynamics of many complex systems, not only in natural sciences but in economical 
and social contexts as well, is usually presented in the form of time series. These series 
are frequently separated by random events which, in spite of their randomness, show some 
*uctu.e and apparent universal featnre. flBJlflQ. During the last few yea. there 
have been endeavors to explain the sort of actions involved in interhuman communication [1|, 
Q]. According to this framework, decisions are taken based on a queuing process and are 
aimed to be valid for a wide range of phenomena such as correspondence among people, the 
consecutive visits of a web portal or even transactions and trading in financial markets ^. 
The main conclusion of these studies is that, in order to reproduce the empirical observations 
as well as to give reason of the heterogeneous nature of outgoing tasks, the timing decision 
has to adopt a rule of non-trivial priority. Otherwise, the implementation of, for instance, 
the simple rule: "first-in-first-out" leads to Poissonian timing between consecutive outgoing 
events and this seems to deviate from many empirical observations. 

One convenient frame to approach these phenomena is provided by the Continuous Time 
Random Walk (CTRW). Within this frame one is basically concerned with the appropriate 
description of the so-called pausing-time density (PTD), which gives the probability 
of having a certain time interval t between two consecutive events. Many empirical PTD's 
present long-tailed profiles suggesting a self-similar hierarchy in the entire probability dis- 
tribution. Following this indication some authors js] claim that the slow decay of the PTD 
obeys a power-law ip(t) ~ whose exponent is almost universal in the sense that it seems 
to adopt only two different values 6 = 1 and 5 = 3/2 [s]. In the next section we will present 
a simple approach which gives a power-law reproducing these exponents. 

Besides the PTD which doubtlessly provides maximal information on interevent statistics, 
the deep structure of the fractal hierarchy is perhaps more easily unveiled by looking at the 
g-moments of the interevent times instead of solely observing the PTD tails. One is thus 
able to answer questions such as whether the process is monofractal or multifractal and 
if there eventually exist different regimes depending on the value of q (the order of the 
g-moment). This information obtained from data can afterwards guide us to find out the 
main ingredients of a more refined theoretical model for human decision dynamics. This is 
certainly the chief motivation of this work. 
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Herein we propose an alternative framework to the existing ones -which are basically 
based on queuing processes- but that it still considers the heterogeneous nature of the 
executed tasks. Within our approach it is possible to deal with analytical expressions, not 
only simulations, and we believe we provide good tools to describe the more subtle structure 
arisen from g- moments. 

The approach we propose has its roots in physics and is reminiscent of Mixture of Distribu- 
tions Hypothesis in Finance that can be traced back to the 1970s Q], the variational principle 
of energy dissipation distributions at different timescales in turbulence in the 1990s |9|, the 
superstatistics and nonextensive entropy [1^. In fact, the PTD ip{t) was first introduced 



within the CTRW model which was originally established by MontroU and Weiss [ill. Il2l. Il3|. 
Under this very general setting, the present development has been inspired by the work of 
Scher and Montroll |l4^ who in 1975 proposed the so-called "valley model" to describe the 
power-law relaxation of photocurrents created in amorphous (glossy) materials. We shall 
use the same idea but in a completely different background. 

The paper is organized as follows. In Sect. [Ill we present the fundamentals of Scher and 
MontroU's model and apply it to explain the emergence of long-tailed distributions in the 
pausing-time statistics. In Sect. Illll we address the question of the moments of the interevent 
times and obtain the conditions for the multifractal behavior of such moments. In Sect. IIVI 
we test multifractality on large financial data sets. Conclusions are drawn in Sect. [Vl and 
some technical details are in the Appendix. 



II. THE VALLEY MODEL AND THE PAUSING-TIME DENSITY 

Scher and MontroU's "valley model" proposes a conditional PTD ip{t\6) as the starting 
distribution. This conditional density accounts for the probability that a given carrier is 
trapped during a time interval t within a potential well of depth e. After this time interval 
has elapsed the carrier jumps to another potential valley. It is next assumed that the energy 
e is a random variable described by a density p{e) 0, Q. We thus have a "superstatistics" 
with the unconditional pausing-time density ip{t) given by 

/oo 
me)pie)de. (1) 
-oo 
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Q 

(2) 



The conditional PTD is assumed to be the simple exponential (Poisson) form 

T{e) [ T{e) 

This choice is quite reasonable since for a given e the emerging statistics is homogeneous 
because all occurrences have the same origin and in consequence they enjoy an identical 
characteristic time scale t{6). Scher and Montroll also assume that the relationship between 
the random energy e and the characteristic time of the distribution is given by the simple 
exponential form: 

T{e) = Toe^' (ro,/5>0), (3) 

where tq = r(0) and (3~^, a fundamental constant of the model, is measured in units of 
energy. We should note that in Scher-MontroU's approach (3~^ = KbT is the thermal 
energy of the environment at temperature T {Kb is the Boltzmann constant). 

We remark at this point that the valley model is consistent with the most basic properties 
of a queuing process recently addressed by Vazquez et al jsl. Indeed, in that process a set 
of incoming messages, or tasks, arrives at random. To these messages a certain priority 
labeled by e is attached. The execution time of a given task with priority e is described 
by the conditional density ijj{t\e). In the most general setting e is also a random variable 
characterized by a density p{e). We are thus faced again with the "superstatistics" mentioned 
above since the timing of the outgoing tasks is governed by the unconditional PTD ilj{t) given 
by Eq. (P). 

In the simplest case of a "first-in-first-out" queue the priority e = fi has the same value 
for all tasks, hence p{e) = 6{6 — p) and the unconditional PTD reads 

t 

r(/i) 

which is a Poissonian density with a single characteristic time scale (the mean time between 
consecutive outgoing events) r(yu) = tqc^^. In terms of decision theory, the situation is 
comparable to that of having no priority protocol at all . 

Another particular situation would be to assign priorities in a uniform random manner 
with p{e) = 1/(2 A) where possible values of e are constrained inside the interval [—A, A]. 
In this case 



^it) = -j-^ exp 
r(/i) 



(4) 
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and from Eqs. 



we can easily see that 



2Apt 



exp — 



T4. 



— exp 



(5) 



where r± = roexp(±A/3) (see Fig. [T]). The characteristic time scale, Tc, of this model (i.e., 
the mean time between outgoing events = {t)) is straightly obtained using Eq. and 
reads 

cosh Ap 



Observe that, when Ti ^ r_ 



2A(3t 



AP 



exp 



-To- 



(6) 



This truncated power-law with exponent 1 is precisely one of the two universal classes 
suggested in Ref. for queuing processes which emerged when the queue had a fixed 
length. In our model, this simple power-law arises when the random variable e is uniformly 
distributed. 
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FIG. 1: (Color online) The pausing-time distribution for different cases as a function of the t/ro. 
Dashed line provides the Poisson case as given by Eq. Dotted line is the PTD Eq. ([5]) when 
p(e) is uniform (tq = 1, A = 3 /3 = 1). Solid line represents the PTD Eq. ^ when p(e) follows 
the Laplace density ^ {a = 2 and (3 = 1. Hence the power law exponent is 3/2, cf. Eq. ([9])). 



The second universal class proposed by Vazquez et al 5(] is given by a power-law with 
exponent 3/2 which appears in their simulated series when the queue is supposed to have an 
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arbitrary length. Let us show that within our approach we can obtain arbitrary power-laws 
of the form ?/'(t) ~ l/t*^ (5 > 1). To this end we shall assume that the random variable e 
follows the Laplace distribution 

P(^) = :^e-l^l/^ (7) 

where a = {\e\) > 0. Plugging into Eq. ([1]), assuming that the conditional PTD is Poissonian 
(cf Eq. ([2])) and using the exponential form of T{e) given by Eq. ([3]), we get 

i'it) = ^ exp [-Pe - {t/To)e-^' - \e\/a\ de. 

The integral appearing in the right hand side of this expression and that runs over the entire 
real line can be splitted in two integrals over the interval (0, oo) 

^(t) = -^[/W(t) + /(-)(t)], 

where 

/(±)(t) = exp [±(3e - (t/ro)e^^" - e/a] de. 
One can easily see by simple transformations of variables that 

and 

n 

where T{a,z) and 7(0, are incomplete gamma functions [16.]. We therefore obtain the 
exact PTD (see Fig. [1] below) 



(f ) 7(1 + t/ro) + (y) r(l - 1/cr/?, t/ro) 



(8) 



The characteristic time scale corresponding to this density is shown to be (recall that 
Tc = (t) and see Eq. ( IT6l) ) 



To 



This characteristic time will exist as long as (3a < 1, i.e, for the "high-temperature phase". 
Thus the process has a completely different dynamics according to whether (3a < 1 or 
(3a > 1. For in the former case the system has a finite average time between consecutive 
outgoing events, while in the latter case, i.e., for the "low-temperature phase", such a time 



ceases to exist, which means that the long term locahsation for the whole observational time 
has sufficient probability. We can say that in this phase event statistics is dominated by 



rare and extreme events 



17 



18|. 



We may therefore assert that when a = I3~^ the system undergoes a "phase transition" . 
Note that for the carrier currents of Scher-Montroll's model such a phase transition will 
occur when the energy fluctuations, a = {\e\), equal the environment's thermal energy, 
= KbT. 

Let us return to Eq. ([8]). For long times t/rg ^ 1 we can write [3] 

7(1 + l/aP,t/To) ~ r(l + 1/aP) - (t/ro)i/'^^e-*/^o 

and 

r(l-l/a/5,t/ro)~(t/ro)-i/'^^e-*/-«. 

Hence, neglecting exponentially small terms, we get 

r(l + l/a/3) /ro\i+i/-/3 
V'W ~ o ^ T » ^o) (9) 



which is a general power-law of the form with exponent greater than 1 (recall that P 
and a are both positive). Therefore, the appearance in the PTD of long tails in the form of 
a nontruncated power-law is exclusively due to the assumption of a nonuniform distribution 
(in the present case the Laplace distribution Eq. ([7j)) for the random variable e. Note 
that in deriving the power law (Q no other power-laws had to be imposed in intermediate 
stages such as the in the conditional PTD, which still maintains the Poissonian form, nor 
in the characteristic time scale r(e) given by Eq. We also observe that, contrary to the 
uniform distribution where e is limited by the fixed length A, there is now no upper or lower 
bound for the random variable e. We finally remark that one of the cases analyzed in j^, a 
power-law with exponent 3/2, comes up from Eq. (Q when a = 2/ (3. 



III. PAUSING-TIME MOMENTS 

The PTD 'il}{t) provides maximal information about the timing between successive events. 
There are, nonetheless, other quantities that can more easily unveil hidden but relevant 
features of the process such as multifractality. This is precisely the case of the pausing-time 
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moments: 

roo 

Jo 

Within our framework based on a superstatistics with some weight function p{e), g- moments 
are written as 

/oo 
{t''\e)p{e)de. (10) 
-oo 

where (f^ls) are the conditional pausing-time moments defined by 



(f^le) = / t'^ilj{t\e)dt. 

JO 

From Eqs. ([2]) and Q we get 

(t'^|e) = r(l + g)roV^ (11) 



and the moments read 



/oo 
e^^^pie)de. (12) 
-oo 



We remark that exponent q may be negative. However, for the Poissonian density ([2]), the 
conditional moment (f^ls) and, hence (f^), will exist as long as g > — 1. Moreover, for the 
moments to exist, p{e) must decay faster than e'^^^ as e ±oo. 

Before proceeding ahead let us briefly recall that the timing process will show multifrac- 
tality if g-moments behave as 

where /(g) is a nonlinear function of g and L is some suitable scale. When /(g) is linear the 



process is termed as monofractal 191] 



We will now explore the possible emergence of multifractality depending on the choice of 
the weight function p{e). The simplest assumption is p{e) = 6{e — p), that is, there is no 
heterogeneity and all the "valleys" have identical depth p. In such a case we readily obtain 
the monofractal behavior 

{t'') = T{l + q)TSl\ (13) 

where 

/ = e^^. (14) 

We have shown in the previous section that the uniform p{e) and the Laplace p{e) both 
result in long-tailed expressions for the PTD ^/'(t). It is thus natural to ask ourselves whether 
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those two long-tailed distributions (cf. Eqs. ([5]) and (j8])) are also responsible for any multi- 
fractal behavior of the g-moments. 

For the uniform distribution p{e) = 1/ (2A) (—A < e < A) we have 

in = ^ - ' (15) 

where t± = Tq exp(±gA/3). Equation (flSi) describes approximately a monofractal situation. 
The genuine monofractal structure with r+ emerges when r_|_ ^ r_ or equivalently if Af3 ^ 1. 
For the Laplace distribution, p{e) = (l/2cr)e^'^'/°', we have 

2a J-oo 

Now g-moments will exist as long a.s f3q < . In such a case we get 

which is again a monofractal. 

Therefore, none of the two weight functions that result in exact long-tailed PTD's produce 
a multifractal structure. To get multifractality we should go one step further in the formalism 
and consider, for instance, a stretched exponential: 



Pi^) = ; i-T exp 

' 2ar(l + l/a) ^ 



fx 



a 



(17) 



(a > and a > 0). Substituting Eq. (ITTI) into Eq. (fT2l) we get after a simple change of 
variables the unconditional moment in the intermediate form 

where / is defined in Eq. f|T4l) and 

/oo 
expi-\yr + qa/3y)dy. (19) 
-oo 

Looking at this integral we see at once that the convergence is assured as long as a > 1. In 
other words the unconditional moment (f^) exists if a > 1 and also if g > —1 (see above). 
We can integrate Eq. ( |T9i) by means of a power series. In effect 

J(g) = 2 e-y" cosHqa(3y)dy = 2 £ ^-^0^ y'^'e-y'dy. 



n=0 



9 



H6I1C6 

r(l + g) - r((2n+l)/a) 

^' ^ = T(i7^^°^ £ — m — ' (2°) 

which is an exact expression for the g-moment. However, due to the term T{{2n + l)/a) 
appearing in the numerator, the convergence of the series in Eq. ( l20l) is quite slow and 
difficult to evaluate from a numerical point of view. Nevertheless, when a = 2 {i.e., for a 
Gaussian e) the integral in Eq. (fT9!) can be done exactly in closed form with the result: 

(t'') = r(l + g)ro^/''e'^W^ 

which, after defining 

clearly shows for the Gaussian e the multifractal behavior of the g-moment: 

=r(l + g)ro^rL^'. (21) 

When a ^ 2 the integral (fTOjl cannot be done exactly and to elucidate any possible 
multifractal behavior of the g-moments we have to resort to approximations. To this end we 
first define in Eq. (fTOjl a new integration variable x by the change of scale y = {PaY^^°'~^^x. 
We have 

/(g) = A^/" exp|-A[|a;r - qx]jdx, (22) 

where 

A = (/5(t)°/("-^). (23) 

Let us suppose that the fiuctuations of the random variable e, represented by a, are larger 
than "the thermal energy of the environment" f3~^, that is, (3a > 1. In such a case the 
dimensionless parameter A is large and we can safely use the saddle-point approximation 
[2^ for the evaluation of /(g). This is done in the Appendix with the result 

/(g) ~ A(g)exp{6|g|°/("-i)}, (24) 

where ^(g) is given in Eq. (\A9\i of the Appendix, and 

b={a-l) {Pa/aT'^"-^^ > 0. (25) 

Substituting Eq. into Eq. (fT8|) we get an approximate analytical expression for the g- 
moment of the interevent time intervals. We write such an expression in a form that clearly 
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enhances its multifractal (MF) character: 

(t") ~ r(l + g)ro^FLl'?l'^'''""'\ (26) 

where, as before, we have introduced two different scales I and L. The former / = e^^ and 
given by Eq. fll4p is related to dissipation because it depends on the dissipative term fi of 
the weight density ffT7|) . On the other hand, the scale 

L = (27) 

providing the multifractal behavior is the responsible for fluctuations since b defined in 
Eq. (125!) depends on a which, in turn, is related to the variance of p{e); the latter given 
by [r(3/a)/r(l/a)]a^. Note that all g-moments considered above obey the normalization 
condition, i.e., they are equal to 1 for g = 0. 

The dissipative and fluctuating scales I and L merge into a single scale when b = 
This equality means that dissipation fi and fluctuation a are linked by 

^ = (28) 

where = (1 — l/a){P /aY^^°'~^\ For a = 2 (the Gaussian case) this relation reads 

/i = (a/2)2 

which is the analog of the usual fluctuation-dissipation relation. This leads us to look at Eq. 
(125]) as the fractional version of the fluctuation-dissipation theorem suitable to the present 
approach. 

Let us finally observe that when the fractional fluctuation-dissipation relation holds the 
monofractal and multifractal parts of the g-moment are both governed by the same scale, 
that is 

(0^r(l + g)ro^L^+l''l"''""^ 



IV. FINANCIAL DATABASE: AN EMPIRICAL ANALYSIS 



We shall now confront our analytical model with empirical data. We focus on moments 
and multifractality and leave for a future presentation extensive tes ting of the PTD's ob 
tained in Sect. [TTlaiid their comparison with previous studies 
29|. 
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TABLE I: Empirical data specifications of the tick by tick intertransaction data used. These 
are futures contracts on German index (DAX), on the Dow Jones American index (DJI), on the 
Pohsh index (WIG20) and on the Foreign Exchange US Dohar-Deutsche Mark (USDM) and US 
Dollar-Euro (EURUS). We also add a single stock: Telefonica (TEF). 



Ticker 


Time Period No. 


of transactions 


DAX 


2007/02/13-2007/06/14 


4 997 027 


TEF 


2006/01/02-2007/08/27 


3010511 


DJI 


2006/03/01-2007/08/27 


3 806 980 


WIG20 


2006/06/19-2007/03/16 


282 007 


USDM 


1993/01/04-1997/07/31 


1 048 590 


EURUS 


2007/08/01-2007/08/27 


4 176 362 



We have decided to apply our approach to financial markets because finance is one of the 
fields where large amounts of data are easily available. In particular we collect tick-by-tick 
data of futures contracts on several indices and also on a single stock (see Table [11). The 
assets chosen have a very diverse nature thus providing wide generality to our analysis. 

Figure [2] displays the empirical g-moments up to order q = 20. The empirical analysis 
shows that all databases adopt a monofractal form for q sufficiently large. To check this we 
define the estimate fj by 

where (tl) is the empirical g-moment of the database labeled by i. We then perform the 
linear regression 

(t'^) 

and obtain that for 10 < g < 20 the estimate Xj is independent of q which proves the 
monofractal character of the g-moments when q is large (see Table [IT] for the specific values 
of Ti). Note that the error in estimation is very small being around 0.03% in all markets. 

Let us now support the monofractal findings for q large and also check the robustness of 
the estimate fj for all data sets. To this end we evaluate the ratio 

^-w- Firing)'. (30) 
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TABLE II: Fitted parameters in seconds of the monofractal case provided by Eq. (j29p in the 
domain W < q < 20. Error in the estimation is very small. 

i dataset Inr^ 

DAX 4.598 ± 0.004 

TEF 5.458 ± 0.004 

DJI 6.144 ±0.005 

WIG20 7.008 ± 0.003 

USDM 5.217 ± 0.003 

EURUS 6.697 ±0.002 




5 10 15 20 

order q 



FIG. 2: (Color online) The normalized moment (t'')/r(l + q) as a function of the order q from six 
empirical financial data sets in a semi-logarithmic scale. Solid lines correspond to regressions with 
functions of the form r'^ for q values between 10 and 20. 

where 6 is an arbitrary parameter identical for all databases while is the parameter 
obtained from the monofractal fit when 10 < g < 20 (see Tab. [TTl)- We have chosen the 
parameter 6 to be the estimated r for the Dow Jones Index, that is, we take 6 = tdji = 
exp(6.144) (in seconds). If the monofractal hypothesis hold, 0j(g) curves would collapse 
into a single straight line (in a semi-logarithmic scale). We show this analysis in Fig. [3] and 
observe that all data sets merge into a single curve for g > 10. Looking more carefully at 
Fig. [3] we also observe a deviation from the monofractal behavior for q < 10. Moreover the 
differences among the data sets become neatly visible for g < 5. 
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order q 

FIG. 3: (Color online) The ratio defined in Eq. ()30p as a function of the order q from six 

empirical financial data sets in a semi-logarithmic scale. Observe the merging of all data sets for 
q > 10. 




EURUS ■ 



2 4 6 8 10 12 14 16 18 
order q 

FIG. 4: (Color online) The function fp^q) defined in Eq. (|3ip as a function of the order q from six 
empirical financial data sets. Observe again the merging of all data sets for q > 10. 



This indeed can be checked in another way than that of Fig. [3] by plotting the quantity 

/F(g) = -r^in((t?)/r(i + g)). (3i) 

q In Ti 

If data were monofractal, the points should be close to 1 independently on the value of q 
but as can be observed in Fig. Hlthis is only true for q > 10. Also observe that the points 
converge to 1 monotonically as the order of the moments increases. Thanks to this, the 
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plot give us some hints on possible multifractal candidates able to fit data for the smallest 
orders of q. The description given by Eq. fl26|) where we assumed a stretched exponential 
for p(e) appears to be a good candidate. Figure O exemplifies these abilities with the stock 
Telefonica and the future contracts on DAX. Table IIIII shows the values of the estimated 
parameters not only for the DAX and Telefonica but also for the rest of datasets. 




01 23456789 10 
order q 



FIG. 5: (Color online) The emphasized plot presented in the semi-logarithmic scale, taken from 
Fig. [21 from six empirical financial data sets for {f^) /T{\ + q). Solid lines provide two different theo- 
retical predictions for the DAX index and Telefonica (TEF) according to multifractal formula ()26p 
with parameters given by Table HIH 

TABLE III: Fitted parameters within the multifractal domain of small g's. Every market has 
different ranges going from g = to g = 3.5. The parameter a appears in general to a smaller 
error involved compared to other parameters. 





q domain 


a 


In(To) ± 


b 


DAX 


0-3.5 


1.85 ±0.13 


-1.5 ±0.4 


0.9 ±0.3 


TEF 


0-3.5 


1.47 ±0.06 


1.45 ±0.12 


0.14 ±0.05 


DJI 


0-3.0 


1.47 ±0.08 


0.34 ±0.16 


0.3 ±0.2 


WIG20 


0-2.0 


1.65 ± 0.06 


1.67 ±0.14 


1.19 ±0.14 


USDM 


0-3.0 


1.89 ± 0.03 


2.82 ±0.02 


0.32 ±0.02 


EURUS 


0-2.5 


2.1 ±0.2 


-4.6 ± 1.1 


3.4 ± 1.1 
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FIG. 6: (Color online) The function /mf(9) defined in Eq. ([33]) as a function of the scaled moment 
order bq^^^"'^^^ from six empirical financial data sets. Observe again the merging of all data sets 
for small q to linear curve with slope equals 1. 

Summarizing we may say that empirical g-moments of financial interevent times clearly 
show multifractal behavior for (approximately) q < 3.5. In this case the analytical model 
expressed by Eq. (l26l) . 

{t'') ^T{l + q)r^l'^L\^\"^'"~'\ (32) 

agrees with empirical data (see Table HTB . On the other hand, for higher values of q the data 
undoubtedly show a clear tendency to monofractality as has been tested in Figs. |3] and HI 
The abilities of the model (|32|) for small q can be checked as well through the function 

/MF(g) = ^ ln((t^)/r(l + q)) - [In(ro) + . (33) 

We can plug into Eq. ( l33il the parameters estimated from each data base and afterwards 
represent /mf(q') as a function of bq^^^°'~^\ In case we observe the merging of all databases 
into a straight line with slope 1 and without independent term, we could thus assert that 
MF model with the stretched exponential for density p{e) is a good candidate (cf. Eq. (|32|) 
and (155]) ). This exercise is done in Fig. O The merging for small q demonstrates the abilities 
of the MF model not only for the Telefonica stock and the futures on DAX but also for the 
rest of financial databases. 

If we want, however, to have all empirical facts in the nutshell of a single formula we 
should generalize Eq. (132]) so as to include the monofractal behavior when q becomes large. 
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FIG. 7: (Color online) The normalized moment (t^)/r(l + q) as a function of the order q from 
Telefonica stock and DAX index in a semi-logarithmic scale. Solid lines correspond to regressions 
assuming the Heuristic Multifractal (HMF) formula (j34p with parameters from Table IIVI while 
dashed lines shows again the Multifractal (MF) model with an stretched exponential superstatis- 
tics (|26p with parameters from Table |TI1 The inset focuses in a shorter domain of q values. 



TABLE IV: Set of the corresponding parameters of the heuristic extension HMF given by Eq. (j34p . 

a In tq + fj,(3 b hi 

DAX 1.91 ±0.03 -3.0 ±0.3 2.5 ± 0.3 0.33 ± 0.02 
TEE 1.78 ±0.02 0.1 ±0.2 1.07 ±0.11 0.20 ± 0.01 
DJI 1.60 ± 0.02 0.18 ± 0.12 0.29 ± 0.05 0.091 ± 0.012 
WIG20 1.96 ±0.05 0.5 ± 0.5 3.3 ± 0.6 0.50 ± 0.05 
USDM 1.69 ± 0.02 2.97 ± 0.06 0.26 ± 0.03 0.115 ± 0.009 
EURUS 2.21 ± 0.03 -9.5 ± 0.7 11.7 ±1.1 0.71 ± 0.03 

The requirements that such a Heuristic Multifractal (HMF) formula has to satisfy are: (i) 
it must obey the normalization condition (i.e., for g = it it should be equal to 1); (ii) 
for small values of q it must reproduce Eq. (!32|) : while (iii) for larger values of q the HMF 
formula must tend to a monofractal form. 
The heuristic formula we propose is: 

(t?) = r(l±g)ro'FL^('?), (34) 
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where 



1 

hi 



l-exp(-6i|g|^/("~^))] (35) 

Note that we have added a fourth parameter, hi > 0, which modifies the scale L = exp(6) by 
a new one exp(6/6i) (cf. with Eq. (126|) ). Equation (IMI) obviously satisfies the normalization 
condition. Moreover, for q small we have 1 — exp feilg]"'^^^"^^-') ~ bi\q\'^^''°'~^'' and we recover 
Eq. ( |32l) . Also for q large 1 — exp (^—bi\q\^^^°'~^^^ ^ 1 and Eq. (IMl) tends to the monofractal 
form: 

(t^) ~r(l + g)ro«/«Ll'?l/^^- 

Figure [7] shows (solid curves) how the HMF formula fits the DAX and Telefonica empirical 
data on the whole range of values of q. This is additionally confirmed for the shorter range 
by the zoom provided by the inset graph. For q < 2 the predictions of MF and HMF 
formulas cannot be distinguished. The value of parameters of the HMF model obtained by 
the fit is given for comparison in Table IIV[ 

The predictions of formula (|34l) have been tested on the available data sets with satisfac- 
tory results by plotting the function 



fnuFiq) = y 



-ln((t^)/r(l + g))-lnro-/i/3 



(36) 



versus feig^'''-""^-'. If the HMF hold, we then would be able to see a merging of all data sets 
along the curve 1 — exp(— x) where x = biq^^^"~^^ (cf. Eqs. f lM|) and fl36l) ). The model also 
allows for studying explicitly the dependence on In L = 6 of transformed g-moment 

"1 



\nm/T{l + q))-\nTo-^I(3 



For doing this we escale the q values in each market i as 



{bi/b\) 



0,;-l 



(ai-l)/(Q-l) 



(37) 



(38) 



where and b\ are these estimated parameters from market i as shown in Table IIV[ The 
unlabeled parameters bi and a concerns a reference market which again corresponds to the 
DJI futures. If the model hold, there would be a linear dependence between the transformed 
g-moments (157|) and InL = 6 across the different markets since the HMF model fl35l) is 
invariant across the markets through 

/hmf(?) = fHMpiQi), 
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FIG. 8: (Color online) The function /hmf(9) defined in Eq. (f36]l as a function of the order biq^^^"'~^^ 
from six empirical financial data sets using parameters of Table llVi We finally observe the merging 
of all data sets for < g < 20 along a curve of the form 1 — exp(— x) where x = biq^^^°^~^\ 
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FIG. 9: (Color online) The transformed g-moment given by Eq. (I37p as a function of InL = 6 when 
q =2 (bottom line), q = 5 and g = 15 (line above) and for six empirical financial data. Solid lines 
verify linear dependence across different markets. 



and where /hmf(q') would be its slope. Figures [HI and [9] show the satisfactory results that 
supports the validity of the Heuristic Multifractal formula. 

We finally mention that the problem of the explicit forms of conditional PTD ^/'(t|e) and 
distribution p{e) which gives the heuristic formula (IMI) according to relation (fTOl) is still a 
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FIG. 10: (Color online) The sojourn probability ^{t) of the Telefonica (TEF) stock. Solid line 
provides the numerical computation of ^{t) for the MF model as given by Eq. (I40p with parameters 
b = 0.12, a = 1.6, and In(ro) + = 1.7. Dashed and dotted lines are respectively the fits with 
the q-exponential (with m = 0.94 it 0.07 and q= 1.3 it 0.2) and the Weibull probabilities (with 
a = 1.53 ± 0.02 and c = 0.459 ± 0.004). 



challenge. We can otherwise check the soundness and self-consistency of our multifractal 
approach by looking at its sojourn probability (i.e, the decumulative probability) 



^(t) 



00 roo 



^{t'\e)p{e)dt'de 



(39) 



and compare it with empirical data. We take the sojourn probability instead of the pausing 
time distribution because verification with empirical data is firmer. Recall that the MF 
model takes the Poisson density ip{t\e) provided by Eq. (j2]) while p{e) obeys a stretched 
exponential as given in Eq. f|T7j) . Substituting these densities into Eq. fl39|) yields 

1 f'^ / , t 



^MF(t) 



exp 



\x\ 



dx. 



(40) 



2r(l + l/a) J~oo - \ ■ ■ roe^ 

The expression needs to be numerically evaluated and for doing this we have taken the 
parameters of Telefonica given in Table IIIII and slightly modify them to improve the fit with 
empirical data. Solid line in Fig. [10] shows the resulting curve and it is there compared with 
the empirical sojourn probability of Telefonica. 



The empirical analysis on the pausing time density and the so iou rn probabi 



data has been extensively studied during the last few years 
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28 



ity in financial 



291 ]. Some recent 
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papers argue that can be described properly by the Tsalhs q-exponential [28|, l29| 

1 

" [l-tm(q-l)t]^/(^-^) 



with q> 1 or the WeibuU distribution 



28 



29| 



^'w(^) = exp {-at") . 

These candidates are also represented in Fig. [10] and the quality of their fits are comparable 
to that of our MF model. A more accurate study among the differences and similarities 
of the MF model (and eventually the HMF) between both the theoretical and empirical 
distributions is certainly necessary. However, we leave a more complete study of the sojourn 
probability for a future work. 

V. SUMMARY AND CONCLUSIONS 

In this work we have extended the original CTRW formalism, within the frame of Scher- 
Montroll's valley model, to furnish an analytical treatment for the statistics of interevent 
times. The model developed has been tested to financial time series, although the analysis 
is applicable to the broader area of interhuman communications. 

The approach presented consists in obtaining the PTD ilj{t) and the g-moments {f^) of 
the interevent time intervals through a random variable e described by a probability density 
p{s). The nature of this hidden variable depends on the problem at hand. In the original 
work of Scher and Montroll e represented the depth of the potential well where carriers were 
trapped. In other contexts, such as queuing processes, e may represent the priority assigned 
to an incoming task and for financial markets we are exploring the possibility that e would 
be related with transaction volumes, market depth or bid-ask spread 

Whatever the case, the overall approach assumes an expression for the conditional PTD 
ip{t\e) governing the timing of incoming events (charged carriers, messages, news, etc.). If 
these incoming events are supposed to arrive at random the natural choice for the conditional 
PTD is the Poisson distribution as given in Eq. ([2]). A second assumption is that for a given 
e the mean time between consecutive events, T{e), depends on the hidden variable e through 
the simple exponential form expressed by Eq. ([3]). Finally, in terms of the probability 
distribution of e the unconditional PTD and the g-moments (which both refer to executed 
tasks or outcoming events) are respectively given by Eqs. ([1]) and (fTOjl . 
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With these simple ingredients we have been able to obtain long-tailed PTD's and multi- 
fractal g-moments. Thus, for instance, for a Laplace density p{e) we have 

m is> 1), 

which agrees with many empirical observations of diverse phenomena from queuing theory 
[5] to finance 27]. 



Regarding moments the choice of a stretched exponential as the probability density for e 
leads to a multifractal behavior of the form 



where L is a conveniently chosen scale. 

We have tested the multi: 
of tick-by-tick data (see 
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ractal behavior of intertransaction times on large financial sets 
, [35I for multifractal analyses in other financial settings). 
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The overall conclusion is that g-moments are multifractal for small values of g (g < 5), while 
for larger orders (f^) becomes monofractal. A more refined but heuristic analytical formula 
has also been proposed which fits the whole range of empirical g-moments. Nevertheless, 
the problem of the explicit forms for both the conditional PTD ip{t\e) and the density p{e) 
resulting in the heuristic expression is still a challenge. 

Let us finish by noting that in some places around the paper we have highlighted some 



thermodynamic similarities in our method. In fact the mu 



;ifractal approach we have herein 



361] . We will develop this analogy 



developed is feasible of a thermodynamic interpretation 19 , 
in a future work. 
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APPENDIX A: APPROXIMATE EVALUATION OF I{q) 

We want to evaluate the integral (12^ which we write in the form 

I{q) = A^/" / e-^^^^'^dx, (Al) 



00 
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where 

h{x) = - qx, (A2) 
a > 1). For A large we can employ the saddle-point approximation or Laplace's method 



201 1 ■ Expanding h{x) = h{xo) + (l/2)/i"(xo)(a; — xq)^ + 0[{x — xq)^] and performing the 



resulting Gaussian integral we obtain 



/(g) = A^/" 



27r 



1/2 



-[A/i(a;o)+0(A-i/2)] 



Xh"{xo)_ 

where Xq is the minimum of h{x). That is, Xq is the solution to 

h'{xo) = a|xo|""^sgn(xo) - g = 0, 



(A3) 



I.e., 



a|xo|° ^sgn(xo) = q. 



(A4) 



But q = |g|sgn(g) and we rewrite Eq. (]A4p as 



a\xo 


a-l 









sgn(g) 



sgn(xo) ' 

Since the right hand side of this equation is positive (recall that a > 0) then necessarily 
sgn(a;o) = sgn(g). Hence 



and the two extremes of h{x) are 



when g > 0, and 







l^^ol = 


(?) 






Xo = 


(g/a) 



l/(a-l) 



(A5) 



xo = -(-g/a)^/(°-^) 



when — 1 < g < 0. 

Taking into account that the second derivative h"{xo) = a{a — l)|a;o|°'~^, i.e. 

(Q-2)/{a-l) 



h"{xo) = a{a-l) I ^ 



> 



(A6) 



is always positive we conclude that Xq is indeed a minimum of h{x). 

On the other hand, recalling that Xq has the same sign as g we can write h{xo) in the 
form 

h{xo) = \xo\°' - \q\\xo\, 
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and using Eq. (lASP we have 



Collecting terms we finally obtain 



h{xo) = -{a - I) r-^] . (A7) 



/(g) ^ A{q) exp , (A8) 

where 



A{q) = A^/" 

and (see Eq. (|23|)) 



27r 



Aa;(a; — 1) 



1/2 /|^|x (2-Q)/2(a-l) 

[a 



(A9) 
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